High dimensional properties of quenched noise growth models. 
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We discuss the behavior of bounded slope quenched noise invasion models in high dimensions. We 
first observe that the roughness of such a steady state interface is generated by the combination of 
the roughness of the invasion process Xc and the roughness of the underlying interface dynamics. 
In high enough dimension we argue that Xc decreases to zero. This defines a critical dimension for 
the problem, over which it reduces to the correlated annealed dynamics, which we show to have 
the same roughness as the annealed equation at five dimensions. We argue that on the Cayley tree 
with one additional height coordinate the associated processes are fractal. The critical behavior is 
anomalous due to strong effects of rare events. Numerical simulations of the model on a Cayley tree 
and high dimensional lattices support those theoretical predictions. 



Kinetic roughening of interfaces during nonequilibrium 
growth has attracted a lot of attraction lately. Initially 
the focus of the theory was on interfaces which are driven 
by random noise which does not depend on the position of 
the interface (annealed noise) A challenging problem 
motivated by experiments is what is the effect of disorder 
in the growth medium (quenched disorder). Examples 
of systems with such disorder are invasion of liquid in 
porous media and magnetic domain wall growth. It was 
discovered experimentally that growth in systems with 
quenched disorder generates exponents which are differ- 
ent significantly from the known annealed exponents. |^] 

In growth problems there is usually a marked distinc- 
tion between the behavior of the quenched noise problem 
and the annealed versions, and they define different uni- 
versality classes. In particular the roughness exponent x 
of the quenched growth is different from Xa, the rough- 
ness of the annealed model. The roughness exponent, 
X, is defined hy W = L^, where W is the system width 
and L is the size of the system. An important example 
is bounded slope growth, which may be described by the 
KPZ equation for the annealed case |]. The quenched 
model §J|,|,| has a different roughness. Growth in 
quenched disorder is characterized by the existence of 
a critical force level. When interfaces are driven at this 
critical force, critical behavior appears, in the sense that 
interfaces become self-afiine, and have infinite correlation 
length, associated with processes of unbounded size. 

In a set of recent papers |||,^ we pointed out that the 
main distinction between growth in annealed disorder 
and growth in quenched randomness is that the growth 
in quenched randomness can be decomposed into consec- 
utive associated processes (APs), which are defined as a 
series of activations below a certain threshold value /o 
of the quenched noise field. These APs are connected 
sets of sites which are delimited by a network of block- 
ing surfaces. The APs have an accompanying roughness 
exponent Xc which is defined by the scaling of the height 



of the APs with respect to their radius. It is usually as- 
sumed that Xc defines the roughness of the interface [^,0 ; 
actually, it is only a lower bound for the overall roughness 
of the surface, not necessarily identical with it 

In this study we concentrate on discrete bounded slope 
models which are believed to belong to the quenched 
KPZ universality class In these models the slope of 

the growing surface is constrained to be below a thresh- 
old value. We discuss the behavior of those models in 
high dimensions. We show that Xc should go to zero in 
high enough dimension. We argue that the APs should 
become fractal objects, with a fractal dimension which re- 
mains bounded in all dimensions. To test this we present 
numerical analysis of systems up to six dimensions and 
on a Cayley tree with a height. We find that there is 
a nontrivial transition between the high dimensional be- 
havior and the Cayley tree. We discuss the relationship 
between our conclusions and a recent analysis Q . 

Analyzing the overall interface growth, the distinction 
between the annealed and the quenched problems be- 
comes less important as soon as the APs become flat (or 
fractal) . The quenched problem then becomes equivalent 
to the correlated annealed noise problem. We analyze the 
relevance of such noise. We argue that the roughness ex- 
ponents of this problem are the same as the annealed 
KPZ exponents at 5 dimensions and above. 

The discrete Buldyrev-Sneppen model |^,^, which be- 
longs to the quenched KPZ universality class, may be 
defined on any space of the form S = Cx Z"*", where £ is 
an arbitrary substrate lattice, and Z+, the set of positive 
integers in the direction of interface growth. C will be 
defined as a hypercubic lattice R'^ in d dimensions and 
as a Cayley tree in the infinite dimensional limit At 
each point s of 5 there is a random pinning force /(s) 
distributed uniformly between and 1. The interface is 
a graph over C. The growth process is as follows: the 
site of minimal / value is activated by increasing the in- 
terface height at the point, followed by growth at nearest 
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neighbor sites, if the local slope exceeds a threshold value. 

The growth of the interface can be represented as a 
series of associated processes, defined as sets of activa- 
tions below a threshold value /o of the quenched field /. 
An activation above /o is concludes an AP and starts a 
new one. The APs of /2 may comprise several /i APs 
if /i < /2. The largest APs are the fc APs, where fc is 
the critical / value, i.e., the largest to be activated in an 
infinite system. A/ = fc~ fo, measures the criticality of 
the APs. The APs are thus random processes character- 
ized by the number of activations s, and the r.m.s radius 
of the projection on C of the process, r|| . 

The following scaling forms are known to be valid for 
d-dimensional APs: 

P{s) = s-^9{s/Arn, Pin) = rp.9||(r||/A/-'^ii), 

(1) 

and the fractal dimension D of the APs is defined by 




The scaling relations 

z^ll = ly/D T|| - 1 = D{t - 1) (3) 

follow immediately from (|l|) and (||). Another scaling 
relation proved in |^ for 1-dimcnsional APs is 

l = i.||(d+l-r||). (4) 

This relation reduces the number of independent static 
exponents to 2. When D > d the APs are rough, with 
roughness exponent Xc — D ^ d. It will be shown below 
that those relations break down for a Cayley tree. 

It is evident that the interface roughness in this system 
cannot be smaller than the value Xc- This is because from 
time to time the surface identifies with a critical surface 
which has this minimal roughness [||. As an example 
one can consider the 1 + 1 dimensional case. Since the 
APs are related to directed percolation their roughness is 
Xc — 0.63 0. The roughness for the annealed dynamics 
is Xa — 0.5, so that Xc dominates Xa- On the other 
hand we should note that the roughness of the annealed 
KPZ equation is a lower bound on the roughness of this 
system. If Xc = 0, the roughness of the system will be 
defined by a process which will be similar to a correlated 
annealed KPZ dynamics. In other words: 

X > max(xQ,Xc)- (5) 

The roughness will be generally defined by an interplay 
between Xa, Xc and the correlations in the system, and 
can actually be larger than both x's. It has been already 
observed that even in the simplest 1 + 1 dimensional case 
the system roughness seems to be higher then Xc Q • For 
higher dimensions the effects should be much stronger 



since Xc decreases. Therefore, we now focus our attention 
on the dependence of Xc on the lattice dimension d. We 
note that the usual relationship x — ^^/^\\ 
is wrong in the general case, this relation applies only 
to chic- Furthermore if x > Xc processes initiated from 
a flat surface such as analyzed in |0,|l^ will be differ- 
ent from the associated processes we defined here in 
scaling properties and in statistics. 

We first present an argument to show that Xc must 
decrease to on the Cayley tree. Consider a large clus- 
ter of radius > r. We want to compare the number of 
activated sites whose distance from a central site is r to 
the number of activated sites with a distance of r + 1. 
Since a critical process is stationary, i.e., the processes 
activated from sites of distance r + 1 are statistically the 
same as those from distance r, the ratio of the number 
of sites on the two shells should be asymptotically 1 for 
a critical cluster. Therefore the number of sites on the 
critical cluster grows more slowly than exponentially, in 
other words, it is fractal, with a fractal dimension Dao- 
Studying the rules of growth indicates that whenever Xc 
is positive, the APs cannot be fractal, since rough clus- 
ters have a typical width r±_ ~ s^"/^ , which through 
avalanches prevents the formation of holes with diameter 
smaller than r_L. For d < dc the APs may be multiply 
connected, but not fractal. APs with Xc = are objects 
with order 1 thickness, which may fractalize. We can 
therefore conclude that Xc is on the Cayley tree. 

We have shown that Xc decreases from 0.63 to when 
d increases from 1 to oo. It is very natural to presume 
that the decrease of Xc as a function of d is monotonic, 
and this is indeed verified by numerical simulations but 
is not critical to the following argumentation. This sce- 
nario for the dependence of Xc on d motivates a defini- 
tion of an upper critical dimension dc where Xc = 0. dc 
marks the transition to fractal APs. We can state that 
dc < Doo < oo. In fact, additional arguments, involv- 
ing dynamic exponents indicate that dc < 6 [Q. Most 
of the preceding arguments on the behaviour of the AP 
as a function of d may be generalized to other growth 
models such as quenched Edwards- Wilkinson growth. 

Fractal APs in this model will be unusual due to rare 
events in which a single site is activated several times, 
creating large, compact avalanches which are exceptional 
in the overall structure of the AP. These events, although 
exponentially rare, can dominate the tail of the distribu- 
tion of the APs, and thus alter the scaling forms (|^). We 
demonstrate this mechanism for the breakdown of the 
scaling form for P(s) on the Cayley tree. Let s{h) be 
the expectation for the size of an AP given that a single 
site is of height h. For h larger then the average height 
s must grow exponentially, 

s{h) ^ /i^ (6) 

with fj, larger or equal to the coordination number of 
the lattice minus 1. /U is an increasing function of /o 



2 



for fo < fc- On the other hand the probabihty that a 
site in an fo AP reaches a height h should decrease like 
f'^ with / larger than, but of the order of fo, and also 
an increasing function of fo- A rough estimate for the 
probability to activate an AP of size s is 



P{s{h) <s< s{h + l)) - /'•; 
combining with M) we get 



P{s) 



log//log/j-l 



(7) 



(8) 



On the Cayley tree, the distribution of the APs is always 
a power law, with a non-universal /o-dependent expo- 
nent. This is in accord with the findings of This 
behavior is in marked contrast to the normal critical be- 
havior (|^). The difference should be attributed to the 
strength of exponentially rare events, which dominate the 
large s tail of the distribution. The main conclusion (|^) 
does not depend on the assumption of localized distribu- 
tion. A more detailed analysis of bounded slope models 
on the Cayley tree will be given elsewhere pl| ]. 

A similar analysis may be carried out for finite dimen- 
sional systems whenever the APs are flat or fractal, so 
that an unusually large height at one site can modify the 
local structure of an AP. In contrast to the Cayley tree, 
avalanches stemming from such rare events are expected 
to grow only as a power of /i, s ~ h'^, so that the large s 
tail of P{s) should decay as a stretched exponential, 



P{s) ~ s '^exp(— as 



(9) 



instead of a power law as in the Cayley tree. The de- 
pendence of a on fo is yet to be investigated. Eq. (||) 
indicates that the mean cluster size does not diverge near 
fc- However, our previous scaling relations (||) and (^) 
and also simulations in high dimensions indicate that 
(s) ^ l/c ^ /o| with 7 > 1. Hence the transition from 
low dimensional behavior to high dimensional behavior 
is anomalous. We conjecture that the scaling form (||) is 
obtained as an infinite dimensional limit of (^. 

We now turn our attention to the scaling of the growing 
interface as a whole. Above dc the APs are flat, and 
growth becomes similar to the annealed process, with 
the source correlation derived from the distribution and 
structure of the flat or fractal APs. In this regime x will 
remain positive and equal to that of annealed KPZ ||l|,^ 

dth(jc, t) = V^/i(x, t) + A(V/i(x, t)f + 77(x, i), (10) 

where h is the height and the random source rj has spatial 
correlations derived from the statistics of the APs. This 
motivates the definition of a second critical dimension 
di, such that Xa — X for d > di. di rather than dc is the 
relevant critical dimension for the overall roughness. 

When d > dcT T] '\s Gaussian and, and its correlation 
function is derived as follows: Since the APs are not 



correlated (for fo = fc, see Q), the noise at two points 
is correlated if and only if the two sites were activated 
in the same AP. Therefore the source correlation may be 
estimated by the probability of such an event, which may 
be estimated as 



(?7(x)?7(y)} - P(r|| > |y - 



x|)|y-x 



D- 



(11) 



We now analyze the effects of those correlations. 
Let h(k,uj) = J d'^xdt/i(x,i)e*-^+*"* with r)(k,uj) de- 
fined similarly. We also define the response function 
G(k,uj)d(k + -k')S{uj +Lu') = {Sh{k,to)/Sr]{k\uj')). The 
scaling form of G is G(k, uj) — k~^g{u!/k^), and the noise 
correlation is {i]{k,Lj)ri{k' ,u;')) = fc"5(k + k')S{uj + uj'), 
where a = D{t — 2). We define the "bare" height field 
as 



/io(k,w) = G{k,uj)ri{k,uj) 



(12) 



If the roughness of ho dominates the annealed white noise 
roughness Xa, the system roughness x must be larger 
than Xa- The bare correlation is 

{ho{k,co)ho{k',uj')) = G{k,uj)G{k',uj') {v{k,uj)r^{k',uj')) 
- k-^''+"g{uj/k'' fS{k + k')S{u; + cj'), 

(13) 

and the simultaneous correlation is given by 



-2z+a 



dujg{uj/k'f 



(14) 



From the correlation in (|T^) we get an effective bare 
roughness 2xo = max(z — a — d, 0). 
The condition for irrelevance becomes 



2x> z — a — d. 



(15) 



15|) is always satisfied for uncorrelated noise (a ~ 0), 
and also reproduces the known condition for irrelevance 
of the noise correlation in 1 dimension, a > —1/2 [Q. 
Numerical simulations give dc ^ Q (see below), and values 
for D and r given in table 1. Inequality (|5|) is satisfied 
for those values, so di < dc- 

We performed numerical simulations of the Buldyrev- 
Sneppen model in order to examine the previous theo- 
retical scenario. The simulations consisted of growth in 
a hypercubic lattice of dimension 1 to 6, and growth on 
a Cayley tree as an effective infinite dimensional lattice. 
The growth on finite-dimensional lattices was simulated 
using the constant current algorithm starting from a fiat 
surface. The search for the global minimum which is 
required in this approach was coded in an efficient man- 
ner, so that about 10^*^ activations were simulated in each 
case. The main computational constraint, however, was 
memory area, which dictated lattices of small sides, down 
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to 16 in 6 dimensions, where there was a small but ob- 
servable scaling range. The simulations were performed 
with periodic boundary conditions. 

The exponents T|| , 7 and D — d were measured. They 
are presented in table 1. These exponents are not in- 
dependent, and 7, for example, is obtainable from the 
other two using relations (0,^). The calculated values of 
7 agree with the measured values of 7, thus providing a 
consistency check on the simulations. (A possible source 
of deviations is a weak size dependence of the exponents.) 

It is evident from table 1 that the APs become flat at 6 
dimensions, so dc ^ 6. The measured Xc becomes equal 
to numerical measurements of Xa fl^ at 4 dimensions 
within numerical accuracy. Since Xc ^ 1 at 5 dimensions, 
it is probably safe to use the criterion ( |l5|) for relevance 
of noise correlation for KPZ roughness. The criterion is 
satisfied so that 4 < < 5. Independent measurements 
of the overall interface roughness indicate that this 
is indeed the case, see table 1. There is no numerical 
indication of an upper critical dimension above which 
the exponents don't change. 

Simulations were performed also of the Cayley tree ver- 
sion of the model as a limit of high dimensional systems. 
Since clusters are expected to be fractal, of a volume 
growing as with the cluster length it is very inef- 
ficient to simulate the process on the full Cayley tree. 
Therefore we simulated the processes on a background 
of a flat blocked surface punctured at one point. As dis- 
cussed above these processes might have different scaling 
than APs on the Cayley Tree because of the strong depen- 
dence of the distributions on height fluctuations [ pl]p4| 
but the general predictions still apply. The clusters were 
grown with a constant driving force, until either the sur- 
face became blocked, or the size of the cluster exceeded 
a certain memory limit (about 2 x 10^ sites). 

To examine the distribution of AP sizes we simulated 
10^ clusters for each value of the driving force /o, on 
a Cayley tree with coordination number 3. As expected, 
we find that the transition is anomalous in the sense that 
the probability distribution function is a power law for 
all values of /o < /c, of the form with no apparent 
cutoff. The value of /c was calculated by estimating the 
number of infinite avalanches initiated by the process. 
An equivalent definition is that t(/c) = 2. It is found 
that fc — .141 ± .001. The measured value of r decreases 
with /o and can be fitted with t(/o) — 2 ^ {fc — foV^^- 
This is the same dependence as in the different Cayley 
tree realization used in 0. The reason for this seems 
to be that the APs are flat so that the differences be- 
tween the models become irrelevant. The relationship 
between the models will be discussed elsewhere [|l^ . The 
measurement of D is more difficult, but a significant de- 
pendence on /o could be observed. The measured value 
at fc was D = 6.5 ± 0.5. This translates to twice the 



ture of quenched noise growth as a function of dimen- 
sionality. We identified the role of the underlying KPZ 
dynamics in determining the roughness of the growing in- 
terface. The concept of AP introduced in Ref remains 
useful in this analysis. The main idea is that those pro- 
cesses possess an intrinsic roughness Xc and that they 
interact via annealed KPZ dynamics. The system rough- 
ness X is a result of this interplay. We show that Xc 
decreases with dimension to at d = dc, and moreover 
the fractal dimension of the APs saturates to a finite 
value -Doo- Hence, above dc the random sequence of APs 
is equivalent to correlated noise acting on annealed KPZ 
equation. Analyzing this effect we show that above 4 di- 
mensions X = Xa- Those general considerations apply to 
any of the various models for quenched disorder discussed 
in the literature ||l|,|ll[. Finally, we have shown that the 
critical behavior of the APs is anomalous on the Cayley 
tree, and made some conjectures as to the transition to 
this behavior from finite dimensions. 
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value in high-dimensional hyper-cubic lattices |13|. 
In this paper we reached an understanding of the struc- 
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TABLE I. Scaling exponents of interface growth, as mea- 
sured from the resuhs of numerical simulations, except 7caic 
which was calculated from Xc and using ^ and (^), and 
X- The numbers in parenthesis are errors in the last digit. 



